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ABSTRACT 

Context. The source XI822-371 is a low-mass X-ray binary system (LMXB) viewed at a high inclination angle. It hosts a neutron 
star with a spin period of ~ 0.59 s, and recently, the spin period derivative was estimated to be (-2.43 ± 0.05) x 10~ 12 s/s. 

Aims. Our aim is to address the origin of the large residuals below 0.8 keV previously observed in the XMM/EPIC-pn spectrum of 
XI822-371. 

Methods. We analyse all available X-ray observations of XI822-371 made with XMM-Newton, Chandra, Suzaku and INTEGRAL 
satellites. The observations were not simultaneous. The Suzaku and INTEGRAL broad band energy coverage allows us to constrain 
the spectral shape of the continuum emission well. We use the model already proposed for this source, consisting of a Comptonised 
component absorbed by interstellar matter and partially absorbed by local neutral matter, and we added a Gaussian feature in absorp¬ 
tion at ~ 0.7 keV. This addition significantly improves the fit and flattens the residuals between 0.6 and 0.8 keV. 

Results. We interpret the Gaussian feature in absorption as a cyclotron resonant scattering feature (CRSF) produced close to the 
neutron star surface and derive the magnetic field strength at the surface of the neutron star, (8.8 ± 0.3) x 10 10 G for a radius of 10 
km. We derive the pulse period in the EPIC-pn data to be 0.5928850(6) s and estimate that the spin period derivative of XI822-371 
is (-2.55 ± 0.03) x 10~ 12 s/s using all available pulse period measurements. Assuming that the intrinsic luminosity of X1822-371 is 
at the Eddington limit and using the values of spin period and spin period derivative of the source, we constrain the neutron star and 
companion star masses. We find the neutron star and the companion star masses to be 1.69 ±0.13 M 0 and 0.46 ±0.02 M 0 , respectively, 
for a neutron star radius of 10 km. 

Conclusions. In a self-consistent scenario in which XI822-371 is spinning-up and accretes at the Eddington limit, we estimate that 
the magnetic field of the neutron star is (8.8 ± 0.3) x 10 10 G for a neutron star radius of 10 km. If our interpretation is correct, the 
Gaussian absorption feature near 0.7 keV is the very first detection of a CRSF below 1 keV in a LMXB. 
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1. Introduction 

Neutron stars (NS) are thought to be born with magnetic fields 
(B-fields) above ~ 10 12 G. Direct measurements of the strength 
of these fields come from the detection of cyclotron resonant 
scattering features (CRSFs). In accreting NS, X-ray pulsations 
and CRSFs are common in systems containing a young high- 
mass companion (HMXBs), with typical strengths of the NS 
B-field between 10 12 —10 13 G, whereas in binary systems con¬ 
taining a low-mass companion (LMXBs) pulsations are detected 
only in a small fraction of systems, and no CRSF has been de¬ 
tected to date. The most likely explanation is that in such sys¬ 
tems, the NS B-field is sufficiently decayed in the course of 
its evolution to a value that, at accretion rates corresponding to 
10 35 —10 38 erg s- 1 , the corresponding magnetospheric radius be¬ 
comes smaller than the NS radius. 

In LMXBs when pulsations are detected, the inferred NS 
B-fields are of the order of 10 8 —10 9 G, which is about three 
orders of magnitude less than the typical values for HMXBs. 
Intermediate values of the NS B-field, between these two ranges, 
are uncommonly observed, most probably for evolutionary rea¬ 
sons. However, some notable exceptions exist, such as (i) the 11 
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Hz pulsar IGR J17480-2446 (iPanitto et al.ll201 ll) . whose NS B- 
field was estimated in the range 2 x 10 s -2.4 x 10 10 G; (ii) the 
2.1 Hz X-ray pulsar GRO J1744-28, with an estimated NS B- 
field of ~ 2. 4 x 10 11 G dCuill 19971) : and, finally, (iii) the 1.7 Hz 
X1822-371 dJonker & van der Klisil^Om ). 


Despite the small sample that it belongs to, the peculiarity of 
XI822-3 71 s till stands out . Analysing the RXTE data of XI822- 
371. Ijonker & van der Klisl (1200 ll ) detected for the first time a 
coherent pulsation at 0.593 s associated with the NS spin period 
and inferred a spin period derivative of (-2.85 + 0.04) x 10 12 
s s- 1 . Analysing RXTE data from 51 976 to 52 883 MJD, 
iJain et all d2010R constrained the spin period derivative better, 
finding (-2.481 ±0.004) x 10~ 12 s s ~*. TTie ephemeris of XI822- 
371 has recently been updated by llaria et all (1201 ll) . who esti¬ 
mated an orbital period of 5.5706124(7) hr and an orbital pe¬ 
riod derivative of (1.51 + 0.08) x 10~ 10 s s _I when analysing 
X-ray data spanning_30 years. A similar sample of data was 
analysed by iBurderi et al.l (1201 Ot) who obtained an orbital pe¬ 
riod deri vative of (_L50 + 0.07) x 10~ 10 s s -1 . In an independent 
paper, IJain et all ( 2010 ) obtained an orbital period derivative of 
(1.3 ± 0.3) x 10~ 10 s s _1 from X-ray d ata a nd, studying the opti¬ 
cal and UV data of X1822-371. lBavless et all ( 2010 ) derived the 
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new optical ephemeris for the source finding an orbital period 
derivative of (2.1 ± 0.2) x 10~ 10 s s _I . 

iBurderi et al.l (1201 Oi) show that the orbital-period derivative 
is three orders of magnitude larger than what is expected from 
conservative mass transfer driven by magnetic braking and/or 
gravitational radiation. They conclude that the mass transfer 
rate from the companion star is between 3.5 and 7.5 times the 
Eddington limit (~ 1.1 x 10 18 g s -1 for a NS mass of 1.4 M Q and 
NS radius of 10 km), suggesting that the mass transfer has to be 
highly non-conservative, with the NS accreting at the Eddington 
limit and the rest of the transferred mass expelled from the sys¬ 
tem by the radiation pressure. iBavless et al.l (120l~0h show that the 
accretion rate onto the NS should be ~ 6.4 X 10~ 8 M© yr~' in a 
conservative mass transfer scenario, again suggesting a highly 
non-conservative mass transfer. 

The large orbital period derivative is a clear clue that the 
intrinsic luminosity of XI822-371 is at the Eddington limit, 
which is almost two orders of magnitude higher than the ob¬ 
served luminosity (i.e. ~ 10 36 e rg s -1 , see e.g . Hellier&Mason 
1989; iHeinz & Nowakl 1200 it iParmar et all 120001: llaria et al 


2001). This is also supporte d by the ratio Lx!L opt of X1822- 


8 x 10 16 G, while for a luminosity of the source of ~ 10 38 erg 


, it assumes a more conceivable value of 8 x 10 10 G. 


found that large residuals are present in the 
EPIC-pn spectrum below 0.9 keV and fitted those residuals by 
adding a black-body component with a temperature of 0.06 keV, 
although they suggest that further investigations were needed to 
understand the physical origin of this component. 

Recently, ISasano et al.l (12014t) have analysed a Suzaku ob¬ 
servation of X1822-371 in the 1-45 keV energy range. The au¬ 
thors detect the NS pulsation in the HXD/PIN instrument at 
0.5924337(1) s and inferred a s pin pe riod derivative of (-2.43 ± 
0.05) x 10~ 12 s/s. ISasano et al.l (i2014l) also suggest the presence 
of a CRSF at 33 keV and inferred from this value a NS B-field 


1 The names of the cited 
sistent with those a dopted 
age XSPEC ( jArnaui 


spectral models are con- 
in the spectral fitting pack- 
3 119961) . The models are described in 


of ~ 2.8 x 10 12 G and a luminosity of the source of ~ 3 x 10 37 
erg s _1 . 

In this work we determine the NS spin period of XI822- 
371 during the XMM-Newton observation and derive a new es¬ 
timation of the spin period derivative, also taking all the mea¬ 
surements of the spin period reported in literature into account, 
including our derived value. We analyse the combined spectra 
of XI822-371 obtained with XMM-N ewton, Chandra (the same 
data sets as analysed b y llaria e t al.l 20131). Suzaku (the same 


data set as analysed by ISasano et al. 


2014 . and INTEGRAL. 


We show the presence of large residuals close to 0.7 keV, while 
we do not find evidence of a cyclotron feature at 33 keV, unlike 
what has been suggested bv ISasano et al.l ( 2014 ). Moreover, we 
show that a CRSF at 33 keV would not be consistent with the 
evidence that the NS in X1822-371 is spinning up. Fitting the 
residuals near 0.7 keV with a CRSF centered at 0.72 keV, we de¬ 
termine a NS B-field (surface) strength between 7.8 x 10 10 and 
9.3 x 10 1() G for a NS radius ranging between 9.5 and 11.5 km. 


2. Observations 


371. Hellier & Mason ( 1989 1) showed that the ratio Lx/L op , for 
XI822-371 is ~ 20, a factor 50 smaller than the typical value of 
1000 for the other LMXBs. This suggests that the intrinsic X-ray 
luminosity is underest i mated by at least a factor of 50. Finally, 
Ijonker & van der Klis ( 2001 :) show that for a luminosity of 10 36 


2.1. The Suzaku observation 


erg s , the NS B-field strength assumes an unlikely value of 


Recently when analysing an XMM-Newton obs ervation of 
XI822-371 and using RGS and EPIC-pn data, llaria et al.l (1201 3l) 
fitted the X-ray spectrum of this source to a model consisting 
of a Comptonised component CompTlQ absorbed by interstellar 
neutral matter and partially absorbed by local neutral matter. The 
authors took the Thomson scattering of the local neutral matter 
into account by adding the cabs 1 component and imposed that 
the equivalent hydrogen column density of the cabs component 
is the same as the local neutral matter. The adopted model is 
similar to the one previously used bv llaria et ah jj2001 kto fit the 
averaged BeppoSAX spectrum of X1822-371 , |laria et akl (1 201 3l) 
suggest that the Comptonised component is produced in the in¬ 
ner regions of the system, which are not directly observable. The 
observed flux is only 1% of the total intrinsic luminosity, the 
fraction scattered along the line of sight by an extended opti¬ 
cally thin corona with an optical depth r - 0.01. This scenario 
explains why the observed luminosity of the source is ~ 10 36 
erg s _1 , while the orbital period derivative suggests an intrinsic 
l uminosity of XI822-371 at the Eddington limit. Furthermore, 
llaria et akl d2013l) 


The X-ray satellite Suzaku observed XI822-371 on 2006 
October 2 with an elapsed time of 88 ks. Both the X-ray Imaging 
Spectrometers (0.2-12 keV, XISs: lKovama_etah[ 2007|)and the 
Hard X-ray Detector (10-600 keV. HXD: lTakahashi et al1l200l 
instruments were used during these observations. There are four 
XIS detectors, numbered as 0 to 3. XIS0, XIS2, and XIS3 all use 
front-illuminated CCDs and have very similar responses, while 
XIS1 uses a back-illuminated CCD. The HXD instrument in¬ 
cludes both positive intrinsic negative (PIN) diodes working be¬ 
tween 10 and 70 keV and the gadolinium silicate (GSO) scin¬ 
tillators working between 30-600 keV. Both the PIN and GSO 
are collimated (non-imaging) instruments. During the observa¬ 
tion, XIS0 and XIS1 worked in 1/4 Window option, while XIS2 
and XIS3 worked in full window. The effective exposure time of 
each XIS CCD is nearly 38 ks, and the HXD/PIN exposure time 
is nearly 31 ks. 

We reprocessed the data using the aepipeline tool pro¬ 
vided by Suzaku FTOOLS version 2(0 and applying the latest 
calibration available as of 2013 November. We then applied the 
publicly available tool aeattcor. s!0 by lohn E. Davis to ob¬ 
tain a new attitude file for each observation. This tool corrects 
the effects of thermal flexing of the Suzaku spacecraft and ob¬ 
tains a more accurate estimate of the spacecraft attitude. For 
our observation, the above attitude correction produces sharper 
point-spread-function (PSF) images. With the new attitude file, 
we updated the XIS event files using the FTOOLS xiscoord 
program. We estimated the pile-up fractions using the publicly 
available tool pileup_estimate. slQ by Michael A. Nowak. 
The pile-up fraction refers to the ratio of events lost via grade 
or energy migration to the events expected in the absence of 
pile-up. The unfiltered pile-up fractions integrated over a circu¬ 
lar region centred on the brightest pixel of the CCD and with a 
radius of 105" are 4.6%, 3.9%, 10.4%, and 9.9% forXISO, XIS1, 
XIS2, and XIS3, respectively. The large pile-up fraction in XIS2 
and XIS3 is due to the two CCDs working in full window dur¬ 
ing the observation. To mitigate the pile-up effects in the spectra 


2 See http://heasarc.gsfc.nasa.gov/docs/suzaku/analysis/suzak 
for more details 

3 http://space.mit.edu/ASC/software/suzaku/aeatt.html 


http: //heasarc. gsfc. nasa. gov/xanadu/xspec/manual/XspecModels Httfipl: //space.mit. edu/ASC/software/suzaku/pile_estimate. si 
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extracted from XIS2 and XIS3, we used annular regions, while 
we adopted circular regions with a radius of 105" to extract the 
spectra from XISO and XIS1. Adopting annulus regions with in¬ 
ner and outer radii of 28" and 105", respectively, the pile-up 
fractions are 4.8% for XIS2 and 4.6% for XIS3. The background 
spectra were extracted using the same regions as were adopted to 
extract the source spectra and centred where the influence of the 
source photons is weak (or absent) in the CCDs. The response 
files of the XIS for each observation were generated using the 
xisrmfgen Suzaku tool, and the corresponding ancillary files 
were extracted using the xisarfgen Suzaku tool, suitable for 
a point-like source. Because the responses of XISO, XIS2, and 
XIS 3 are on the whole very similar, we combined their spectra 
and responses using the script addascaspec. The XIS spectra 
were rebinned to have 1024 energy channels. 

We also extracted the PIN spectra using the Suzaku tool 
hxdpinxbpi. The non X-ray and cosmic X-ray backgrounds 
were taken into account. The non X-ray background (NXB) was 
calculated from the background event files distributed by the 
HXD team. The cosm ic X-ray background (CXB) is from the 
model bv lBoldtl (1 1987h . The response files provided by the HXD 
team were used. The GSO data were not used, considering the 
low signal-to-noise ratio above 40 keV. 

2.2. The XMM-Newton observation 

The region of the sky containing XI822-371 was observed by 
XMM-Newton between 2001 March 07 13:12:48 UT and March 
08 03:32:53 UT (Obs. ID. 0111230101)for a duration of 53.8 ks. 
The European Photon Imaging Camera (EPIC) on-board XMM- 
Newton consists of three co-aligned high-throughput X-ray tele¬ 
scopes. Imaging charge-coupled-device (CCD) detectors were 
placed in the focus of each telescope. Two of the CCD detec - 
tors are Me tal Oxide Semi-conductor (MOS) CCD arrays (see 
iTurner et al.ll200ll) . while the third came ra uses pn CCDs (here¬ 
after EPIC-pn, see IStriider et al.l 1200 ll) . Behind the two tele¬ 
scopes that have the MOS cameras in the focus, about half of 
the X-ray light is utilised by the reflection grating spectrometers 
(RGSs). Each RGS consists of an array of reflection gratings that 
diffracts the X-rays to an array of dedicated C CD detectors (see 
iBrinkman et al.lfl998^ iden Herder et al.ll200lh . 

During the observation, MOS 1 and MOS2 camera were op¬ 
erated in fast uncompressed mode and small window mode, re¬ 
spectively. The EPIC-pn camera was operated in timing mode 
with a medium filter during the observation. The faster CCD 
readout results in a much higher count rate capability of 800 
cts/s before charge pile-up become a serious problem for point¬ 
like sources. The EPIC-pn count rate of the source was around 
55 cts/s, thereby avoiding telemetry and pile-up problems. 

Although the RGS and EPIC -pn data products were ex¬ 
tracted and analysed bv llaria et al.l (1 201 3 i). we extracted the data 
products of the RGS and EPIC-pn camera again using the very 
recent science analysis software (SAS) version 13.5.0 and the 
calibration files available on 2013 Dec. 17. We used the SAS 
tools rgsproc, emproc, and epproc to obtain the RGS, MOS, 
and EPIC-pn data products. 

Since the EPIC-pn was operated in timing mode during the 
observation, we extracted the EPIC-pn image of RAWX vs. 
PI to select appropriately the source and background region. 
The source spectrum is selected from a box region centred on 
RAWX=38 with a width of 18 columns. The background spec¬ 
trum was selected from a box region centred on RAWX=5 with 
a width of two columns. We extracted only single and double 
events (patterns 0 to 4) for the source and background spectra 


and applied the SAS tool backscale to calculate the different 
areas of the source and background regions. 

We extracted the MOS 1 source spectrum, adopting a box re¬ 
gion centred on RAWX=317 having a width of 50 pixels. The 
MOS 1 background spectrum is extracted from a source-free re¬ 
gion selected in one of the outer CCDs that collect photons in 
imaging mode during the observation. We extracted the MOS1 
spectrum and the corresponding redistribution matrix and an¬ 
cillary files using the standard recipe^. We also extracted the 
source+background light curve, observing that the average count 
rate during the observation is 15 c/s. The MOS2 source spectrum 
was extracted from a circular region centred on the pixel show¬ 
ing the largest number of photons; the radius of the region is 
640 pixels. A circle with radius of 640 pixels was placed in a 
source-free region to extract the background spectrum. We ex¬ 
tracted the MOS2 spectrum and the corresponding redistribution 
matrix and ancillary files using the standard recipe 0 We also 
extracted the background-subtracted light curve using the SAS 
tool Epiclccorr, and the average count rate is 15 c/s. Since the 
count-rate limit for avoiding pile-up for the MOS cameras is 100 
c/s and 5 c/s for a point-like source in timing uncompressed and 
small window modefl respectively, we expect that the pile-up 
effects are present in the MOS2 spectrum. The RGS1, RGS2, 
MOS1, MOS2, and EPIC-pn spectra have an exposure time of 
53, 51,51,51, and 51 ks, respectively. 


The Chandra satellite observed X1822-371 from 2008 May 20 
22:53:00 to 2008 May 21 17:06:51 UT (Obs. ID. 9076) and from 
2008 May 23 13:20:56 to 2008 May 24 12:41:54 UT (Obs. ID. 
9858) using the HETGS for total observation times of 66 and 
84 ks, respectively. Both observations were performed in timed 
faint mode. The two Chandra observations had already been 
analysed by llaria et akl (120131) . In this work we only used the 
first-order MEG spectrum fsee llaria et alJ 120131 for details on 
the data extraction) because the first-order HEG spectrum has a 
much lower effective area with respect to MEG below 1 ke V, and 
its analysis is not useful for investigating the presence of features 
in the spectrum of X1822-371 at those energies. The first-order 
MEG spectrum has a total exposure time of 142 ks. 


The INTErnational Gamm a-Ray Astrophysics Laboratory 
(INTEGRAL IWinkler et al.l l2003h has repeatedly observed the 
X1822-371 region. We searched the whole IBIS (lUbertini et al.l 
l2003i ) and JEM-X ( Lund et al .l’ 2003 ) public catalogues, select¬ 
ing only pointings (science windows, SCW) with sources within 
six degrees of the centre of the field of view and with exposures 
longer than 500 s to reduce the calibration uncertainties of the 
IBIS/ISGRI dLebrun et al.l l2003i) spectral response. The avail¬ 
able IBIS data set covers the period starting from 2003 March 
21 until 2013 March 21 for a total usable on-source time of 
1 271 ks and an effective dead-time corrected exposure of 874 ks. 
Because of the smaller field of view, the total exposure of JEM- 
XI (camera 1) is 330 ks, while the dead-time-corrected exposure 
is 283 ks. The INTEGRAL data analysis uses standard proce¬ 
dures within the offlin e scie nce anal ysis s oftware (OSA10.0) dis¬ 
tributed by the ISDC ( ICourvoisier et al.ll2003 i. In the catalogue 

5 see http://xmm.esac.esa.int/ sas/ current/documentation/threE 

6 see http://xmm.esac.esa. int/ sas/current/documentation/threE 

7 see http://xmm.esac.esa.int/external/xmm_user_support/docun 


2.4. The INTEGRAL observations 


2.3. The Chandra observations 
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used for the extraction of the IBIS spectra, we have included 
all the sources significantly detected in the total image obtained 
by mosaicking the individual pointings. We exploited a custom 
spectral binning optimised in the energy range 20-100 keV, and 
in the detection of spectral features around 30 keV, weighted the 
time-evolving response function according to the available data 
and excluded the data below 21 keV due to the evolving de¬ 
tector’s low threshold. For JEM-X, we adopted the standard 16 
bins spectrum provided by the analysis software and excluded 
the data below 5 keV and above 22 keV, which are affected by 
calibration uncertainties. @ 

3. Search for the spin period in the XMM data 

We used the EPIC-pn events to search for the spin period. We 
applied the barycentric correction with respect to the source 
coordinates, given by llaria et ail (1201 ll) . using the SAS tool 
barycen; subsequently, we corrected the data for the orbital 
motion of the binary system using the recent X-ray ephemeris 
of X1822-371 derived by llaria et akl (1201 ll) (see Eq. 2 in 
t hat w ork) and the a sin i value of 1.006 lt-s (see Table 1 in 
Llonker & van der Klisll200fh . We selected the EPIC-pn events in 
the 2-5.4 keV energy range and explored the period window be¬ 
tween 0.592384 and 0.593384 s using the FTOOL efsearch in 
the XRONOS package. We adopted eight phase bins per period 
(a bin time close to 0.074 s) for the trial folded light curves and 
a resolution of the period search of 1 x 10~ 6 s. We observed a 
X 2 peak of 41.89 at 0.592884 s, as shown in Fig. Q] We fitted 
the peak with a Gaussian function, assumed the centroid of the 
Gaussian as the best estimation of the spin period, and associ¬ 
ated the error derived from the best fit. We find that the spin 
period during the XMM observation is 0.5928850(6) s, and the 
associated error is at the 68% confidence level. 

Considering that we have seven degrees of freedom, the 
probability of obtaining a^ 2 value greater than or equal to 41.89 
by chance is 5.47 x 10~ 7 for a single trial. In our search we 
adopted 10 3 trials (we span 10 3 s with a resolution of the pe¬ 
riod search of 1 x 10~ 6 ), consequently we expect a number of 
- 5.5 X 10 4 periods with a _y 2 -value greater than or equal 
to 41.89. This implies that our detection is significant at the 
99.945% confidence level. 

Then, we folded the 2-5.4 keV EPIC-pn light curve us¬ 
ing the spin period of 0.5928850(6) s and adopting 16 phase 
bins per period. We used the arbitrary value of 51975.85 MJD 
as epoch of reference. The folded light curve is almost sinu¬ 
soidal (see Fig. 0. Fitting the folded light curve with a constant 
plus a sinusoidal function with period fixed at one, we obtain 
a;y 2 (d.o.f.) of 13.9(13). Since the constant is 23.92(2) c/s, the 
background count rate is close to 1.3 c/s and the amplitude of 
the sinusoidal function 0.17(3) c/s, we estimate that the pulse 
fraction is 0.75 + 0.13% c ompati bl e w ithin 3 <t t o the value of 
0.25 ± 0.06% reported by llonker & van der Klisl (120011) using 
RXTE/PCA data in the same energy band. 

We report in Table Q] the 13 values of the spin period of 
XI822-371 and the corresponding errors previously estimated, 
together with the value found in the present work. The corre¬ 
sponding times are the mean values between the start and stop 
time of the observations in which the spin period was detected, 
the associated errors are one half of the duration of the corre¬ 
sponding observation. After deriving the geometric mean of the 

x The JEM-X2 unit was active only during a limited part of the mis¬ 
sion, so the exposure time is not enough to provide a significant spectral 
constraint. 



Period (s) - Offset 


Fig. 1. Folding search for periodicities in the 2-5.4 keV EPIC- 
pn light curve. We adopt 8 phase bins per period for the trial- 
folded light curves and a resolution of period search of 1 x 1CF 6 s. 
The peak of ^ 2 is detected at 0.592884 s. The horizontal dashed 
line indicate the x 2 value of 18.47 at which we have the 99% 
confidence level for a single trial. 



Spin Phase 

Fig. 2. EPIC-pn folded light curves ofX1822-371 using the fold¬ 
ing period 0.5928850(6) s. The folded light curve is obtained 
using 16 phase bins per period. 

times in Table Q] (column 1) obtaining T mi ,j - 52,257.78 MJD, 
we fitted the spin periods with respect to the times with T mi( i sub¬ 
tracted using ajinear function to estimate the spin period deriva¬ 
tive. Unlike lJain et al.l d 201(j ). we did not obtain a good fit, since 
the;y 2 (d.o.f.) was close to 10 5 (12); however, we obtained a very 
high value of -0.9992 of the Pearson correlation coefficient. The 
reason for the high x 2 value is not clear to us, but it could be due 
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Table 1. Times and corresponding spin periods 


Times 

(MJD) 

Spin period 

(s) 

Ref. 

50,352.9(6) 

0.59325(2) 

1 

50,993.4(6) 

0.59308615(5) 

1 

51,975.9(3) 

0.5928850(6) 

2 

51,976.04(6) 

0.59290132(11) 

3 

52,094.80(6) 

0.59286109(8) 

3 

52,095.73(6) 

0.59286421(12) 

3 

52,432.62(18) 

0.5927922(13) 

3 

52,489.7(9) 

0.5927790(6) 

3 

52,491.61(15) 

0.5927795(11) 

3 

52,503.45(12) 

0.5927737(10) 

3 

52,519.41(18) 

0.5927721(8) 

3 

52,882.15(9) 

0.5926793(15) 

3 

52,883.15(9) 

0.5926852(21) 

3 

54,010.0(6) 

0.5924337(10) 

4 


Note — Referen ces: 1 Jonker & van der Klisl d200 if ). 2 this work, 3 
Jain et al .. ( 20lQh . and 4 lSasano^t all d2014T) . 

to an underestimation of the errors (e.g. small differences in the 
orbital ephemeris used to correct the data) or model complica¬ 
tions (e.g. caused by small fluctuations around an average linear 
trend), or other issues. The detailed investigation of this aspect 
goes beyond the aim of this paper. 

Considering the high value of the x 2 for the linear fit and 
to estimate the error associated to the spin period derivative, we 
fitted the 14 points without the estimated errors and attributed 
the post-fit errors to the best-fit parameters under the assump¬ 
tion that the model is reliable. In this way we at least get an 
estimation of the averaged linear trend of the measured spin pe¬ 
riod with respect to time. Fitting the data again with a linear 
function, we obtain a = 0.592826(6) s and b - -2.20(3) x ICC 7 
s/d, with the errors at 68% confidence level. This implies that the 
spin period derivative, P s , is -2.55(3) X 1CT 12 s/s; this is compat¬ 
ible within three sigmas with the previously reported values_of 
-2.481(4) x 10 12 s/s and -2.43(5)x ICC 12 s/s given bv lJain et~aj] 
d20l0l) and lSasano et al.l (1201 4) . respectively. We show in Fig. [3] 
(top panel) the 14 points and the corresponding linear best fit. 
The corresponding residuals are shown in Fig. [3] (bottom panel). 

Furthermore, we search for the same periodicity in the 
MOS1 data (taken in timing mode); unfortunately, the lower 
statistics with respect to the EPIC-pn data do not allow us to 
detect the periodicity in this data set. 

We note that the spin period obtained from the EPI C-pn 
events is inconsistent with the value reported by I Jain et all 
(120101) using RXTE/PCA observations simultaneous to the 
XMM observation that we analyse in this work (see Tableland 
Fig.0. To confirm the robustness of our results, we reanalysed 
the simultaneuous RXTE/PCA observations (P50048-01-01-00, 
50048-01-01-01, P50048-01-01-02, P50048-01-01-03, P50048- 
01-01-04, P50048-01-01-05, P50048-01-01-06) spanning 2001 
March 7 10:29:26 UT to 2001 March 8 4:47:44 UT. We applied 
the barycentric correction with respect to the source coordinates 
and corrected the data for the orbital motion of the binary system 
as done for the EPIC-pn data. No energy selection was applied to 
the RXTE/PCA events file. We explored the period window be¬ 
tween 0.592384 and 0.593384 s using the FTOOL efsearch in 
the XRONOS package. We adopted eight phase bins per period 
for the trial-folded light curves and a resolution of the period 
search of 1 x ICC 6 s. We observe ax 2 peak of 150 at 0.592884 s, 
as shown in Fig. [4] We fitted the peak with a Gaussian function. 



-2000 -1000 0 1000 2000 
Time-T m|d (d) 


Fig. 3. Top panel: spin period values shown in Table Q] vs. 
time in units of days (see the text). The linear best fit is also 
plotted. The black squares, open green squares, blue diamond, 
and red triangle indicate the spin perio d value s reported by 
iJonker & van der Klisl (1200 lh . iJain et all (1201 Oh . ISasano et al.l 
(120141) . and this work, respectively. Bottom panel: the corre¬ 
sponding residuals in units of 10 6 s. 


assumed the centroid of the Gaussian as the best estimation of 
the spin period, and associated the error derived from the best 
fit. We find that the spin period during the RXTE/PCA obser¬ 
vations is 0.5928846(3) s and that the associated error is at the 
68% confidence level. This result confirms our detection in the 



Fig. 4. Folding search for periodicities in the RXTE/PCA obser¬ 
vations. The peak of^- 2 is detected at 0.592884 s. 
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XMM/Epic-pn data. Finally, we note that the value of the spin 
period derivative obtained above does not significantly change 
when using our value instead of the one reported in Table Q] 


4. Energy range selection for the XMM-Newton, 
Suzaku, Chandra, and INTEGRAL spectra 

We rebinned the MOS1 and MOS2 spectra using the SAS tool 
specgroup to have at least 25 counts per energy channel and 
with an over-sample factor of 5. To verify that the spectra are 
similar and see how the pile-up effects influence the MOS2 
spectrum, we fitted the two MOS spectra in the 0.6-10 keV en¬ 
ergy range using XSPEC (version 12.8.1). We adopted a very 
simple model consisting of phabsj^pcfabs*CompTT 1 , similar 
to the model used by llaria et akl (1200 ll) to fit the BeppoSAX 
broad-band spectrum of XI822-371. The phabs component 
takes the photoelectric absorption by the interstellar neutral mat¬ 
ter into account, it is a multiplicative component defined as 
M(E) = exp[-/V w * cr(£,’)], where cr{E) is the photo-electric 
cross-section (not including Thomson scattering). The free pa¬ 
rameter, Nh , is the equivalent hydrogen column density in units 
of 10 22 atoms cm~ 2 . The pcfabs component takes the pho¬ 
toelectric absorption into account owing to the neutral matter 
near the source, and it is a multiplicative component defined as 
M(E) - f * exp [~Nh pc * cr(£)] + (1 - /), where Nh pc is the equiv¬ 
alent hydrogen column in units of 10 22 atoms cm~ 2 of the local 
neutral matter, and / is a dimensionless free parameter ranging 
between 0 and 1 that takes the fraction of emitting region oc¬ 
culted by the local neutral matter int o acc ount. 

The CompTT component (iTitarchukl 119941) is a 
Comptonisation model of soft photons in a hot plasma. For this 
component, the soft photon input spectrum is a Wien law with a 
seed-photon temperature, kTo, that is a free parameter. The other 
free parameters of the component are the plasma temperature, 
kT e , the plasma optical depth r, and the normalisation N CompTT- 
We used a slab geometry for the Comptonising cloud. We fitted 
the 0.6-10 keV MOS1 and MOS2 spectra simultaneously with 
the aim of estimating the pile-up effects in the MOS2 spectrum. 
We obtained a large ^(d.o.f.) of 1716(540). The two spectra 
are consistent with each other between 0.6 and 7 keV. Above 
7 keV the pile-up distortion is evident in the MOS2 spectrum. 
Furthermore, the presence of two emission lines in the residuals 
at 6.4 and 6.97 keV is evident; finally, both the spectra are not 
well fitted between 0.6 and 1 keV and show large residuals. We 
show the MOS1 and MOS2 residuals in Fig. 0 We sum the 
MOS 1 and MOS2 spectra using a new recipqj. In the following 
the summed spectrum is called MOS 12 spectrum. The MOS 12 
spectrum is rebinned with an over-sample factor of 5. 

We added the first-order spectrum of RGS1 and RGS2 to¬ 
gether using the SAS tool rgscombine; hereafter, the summed 
spectrum is called RGS12. We rebinned the RGS12 to have at 
least 200 counts per energy channel. In the following, we anal¬ 
yse the RGS12 spectrum in the 0.35-2 keV energy range. 

The EPIC-pn spectrum is rebinned using the SAS tool 
specgroup imposing at least 25 counts per energy channel 
and an over-sample factor of 5. To check the consistency of 
the EPIC-pn and MOS 12 spectra, we fit them simultaneously 
adopting the same model described above. We obtain a large 
^f 2 (d.o.f.) value of 3228(466). We show the two spectra (MOS 12 
and EPIC-pn spectra) and the corresponding ratio (data/model) 
in Fig. [6] We observe a large absorption feature in both spectra 



Energy (keV) 

Fig. 5. MOS 1 (black) and MOS2 (red) residuals with respect to 
the model described in the text. 



Energy (keV) 

Fig. 6. MOS 12 (0.6-7 keV; black) and EPIC-pn (0.6-10 keV; red) 
spectra (top panel) and ratio (data/model; bottom panel) with 
respect to the model described in the text. 


at 0.7 keV. A mismatch between the two spectra is evident be¬ 
low 1.5 keV; the EPIC-pn residuals show an instrumental feature 
at 2.2 keV owing to the neutral gold M-edge. Finally, we detect 
the presence of strong emission lines in the Fe-K region. The 
causes of the mismatch between the EPIC-pn and MOS 12 are 
not clear. However, we note that the residuals close to 0.7 keV 
are present in both the EPIC-pn and MOS 12 spectra (see Fig. 
[6] lower panel). We adopt the 0.6-10 keV energy range for the 
seehttp://xmm.esac.esa.int/sas/current/documentation/ilfiPBadphEgiBClmtingiftgfaH;lialinstrumental feature at 2.2 keV in 
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phabs comp onent). We used the abundances provided by 
[Asplund et al. ( 2009 ) and the photoelectric cross-section given 
by Verner et al. ([ 1996b . We took the Thomson scattering of the 


local neutral matter into account by adding the cabs component 
and imposed that the equivalent hydrogen column density of 
the cabs component is the same as the local neutral matter. 
The constant f gives the percentage of emitting region occulted 
by the local neutral matter. Finally, LN and Ed in the model 
indicate all the Gaussian components added to the model to 
fit the several emission lines observed in the spectrum and the 
added absorption edges, respectively. 

Since the XMM-Newton, Suzaku, Chandra, and 
INTEGRAL observations are not simultaneous we left the 
values of the electron temperature and of the optical depth of 
the CorapTT component free to vary independently. The depths 
of the absorption edges added above 7 keV were free to vary 
independently for the XMM-Newton, Suzaku, and INTEGRAL 
spectra, whilst they were tied in the Chandra spectrum to the 
values of the XMM-Newton spectrum because the Chandra 
spectrum extends up to 7 keV. 


Energy (keV I 

Fig. 7. HXD/PIN source spectrum (black) and NXB+CXB spec¬ 
trum (red) are shown. The NXB+CXB spectrum overwhelms 
the source spectrum at energies larger than 36 keV, the energy 
threshold is indicated with a dashed vertical line. 


the EPIC-pn spectrum, we use an absorption Gaussian line with 
the centroid and the width fixed at 2.3 and 0 keV, respectively. 
We fit MOS12 spectrum using the 0.6-7 keV or the 1.5-7 keV 
energy band. 

For the Suzaku/XIS data, we adopt a 0.5-10 keV energy 
range for the XIS0+XIS2+XIS3 (hereafter XIS023) and XIS1 
spectra. We exclude the energy interval between 1.7 and 2.4 keV 
in the XIS023 and XIS1, because of systematic features associ¬ 
ated with neutral silicon and neutral gold edges. We group the 
HXD/PIN spectrum to have 25 photons per channel and use the 
energy range between 15 and 36 keV. We show in Fig. [7] the 
HXD/PIN source spectrum and the summed HXD/PIN NXB and 
CXB spectra (hereafter NXB+CXB spectrum). The NXB+CXB 
spectrum dominates the source spectrum at energies higher than 
36 keV. 

We analyse the Chandra/MEG in the 0.5-7 keV energy range. 
Finally, we analyse the JEM-X and IBIS spectrum in the 5-22 
and 21-60 keV energy band, respectively. The IBIS spectrum is 
background-dominated above 60 keV. 



5. Spectral analysis 

We simultaneously fitted the XMM-Newton, Suzaku, Chandra, 
and INTEGRAL spectra. We add ed a systematic error of 1% 
to take into account that the observations are not simultaneous. 

Initially, we fitted the spectra using the 0.6-7 keV energy range 
for the MOS12. 

W e fitted the spectra using XSPEC version 12.8.1 (see 
lArnaudl 11996 ). Initially, to fit the _continuum emission, 
we adopted the model used by llaria et ail d2013l) . It is 
Ed*phabs*(f*cabs*phabs*(LN+CompTT) + (l-f)* (LN+CompTTD’O jtxvi. The emission lines are fitted with Gaussian components, 
which is a Comptonised component (CompTT in XSPEC) ab- We fixed the energies and widths of the emission lines below 
sorbed by neutral interstellar matter (the first phabs component) 6 keV at their best-fit values because their analysis is not the 
and partially absorbed by local neutral matter (the second aim of this work; a detailed analysis of these lines is reported 


Energy (keV) 

Fig. 8. Residuals with respect to the best-fit models shown in 
Tables E] and |3] The RGS12, EPIC-pn, XIS023, XIS1, MOS12, 
HXD/PIN, MEG, JEM-X, and IBIS spectra are shown in black, 
red, green, blue, magenta, light-blue, yellow, grey, and orange, 
respectively. The data are graphically rebinned. From top to bot¬ 
tom, the residuals with respect to the continuum consist of: 1) a 
Comptt partially absorbed by local neutral matter (large resid¬ 
uals are evident at 0.7 keV); 2) gabs*Comptt with the energy 
of the gabs component close to 0.73 keV; 3) cyclabs*Comptt. 
The MOS12 spectrum covers the 0.6-7 keV energy band. 

Several emission lines are detected and identified with 
N vii, O vii intercombination line, O vm, Ne ix intercombination 
line, Ne x. Mg xi intercombination line. Mg xn, Si xiv, Fe i, and 
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Table 2. Best-fit values of the continuum emission 




MOS12 spectrum in 0.6-7 keV 



MOS12 spectrum in 1.5-7 keV 



XMM 

Suzaku Chandra/MEG 

INTEGRAL 

XMM 

Suzaku Chandra/MEG 

INTEGRAL 

Parameters 

N h (10 22 cm^ 2 ) 


0.126 ±0.003 

Model: 

Comptt 

0.121 ±0.003 


N V (10 22 cm- 2 ) 


4.48 ± 0.08 



4.57 ± 0.08 


/ 


0.633 ± 0.005 



0.638 ± 0.005 


Egdgel (keV) 


7 20 +0-02 

' - - 0.03 



7 20 +0 02 
' • - 0.03 


^Edgel 

0.076 ±0.013 

0.130 ±0.011 0.076 ±0.013 

< 0.084 

0.072 ±0.013 

0.129 ±0.011 0.072 ±0.013 

< 0.084 

EEdge 2 (keV) 


8.42 ± 0.05 

0 13+0.08 
u - -0.10 


8.43 ± 0.05 

r\ i 9+O.O8 

U ' A -0.10 

TEdge2 

0.09 ± 0.02 

0.084 ± 0.013 0.09 ± 0.02 

0.10 ±0.02 

0.084 ±0.013 0.10 ±0.02 

kT 0 (keV) 


0.059 ± 0.006 



0.058 ± 0.006 


kT e (keV) 

3.55 ±0.11 

4.32 ± 0.05 3.2 + !!^ 

4.89 ± 0.08 

3.63 ±0.12 

4.33 ± 0.05 3.3+p | 

9.04 ± 0.09 9.6)° : J 

4.89 ± 0.08 

T 

9.60 ±0.15 

9.09 ± 0.09 9.7+° ! 

6.0 ± 0.2 

9.43 ±0.15 

6.0 ±0.2 

NcompJT (10 2 ) 


6 6 +0 - 3 
° -0.2 



6.5 ±0.3 


Xl/d.o.f.) 


1.21(2496) 



1.10(2466) 





Model: gabs*Comptt 



N H (10 22 cm- 2 ) 


0.100 ±0.006 



0.100 ±0.006 


N Hpc (10 22 cm- 2 ) 


4.54 ± 0.08 



4.61 ± 0.08 


/ 


0.645 ± 0.005 



0.647 ± 0.005 


E Edgel (keV) 


7 20 +0 02 
- 0.04 



7 20 +0 02 
' • - 0.04 


^Edgel 

0.075 ± 0.013 

0.125 ± 0.011 0.075 ±0.013 

< 0.084 

0.071 ±0.013 

0.126 ±0.011 0.071 ±0.013 

< 0.084 

EEdge2 (keV) 


8.43 ± 0.05 

0 i2 +0 - 08 
U ' 1Z, -0.10 


8.43 ± 0.05 

Q J 2 + 008 
U ' 1Z -0.10 

TEdge2 

0.10 ±0.02 

0.081 ±0.013 0.10 ±0.02 

0.10 ±0.02 

0.081 ±0.013 0.10 ±0.02 

Egabs (keV) 


0.73 ± 0.03 


0.72 ± 0.03 

O -,gabs (keV) 


0.14 ±0.03 



0.13 ±0.03 


Tgabs 


0 . 038 +qoo 9 



0 031 +0 013 
u.udiq 009 


kT 0 (keV) 


0.056 ± 0.010 



0.052 ±0.011 


kTf, (keV) 

3 61 +012 

j.ui 010 

4.33 ± 0.05 3.2+° | 

9.00 ±0.10 9.6+° : ^ 

4.89 ± 0.08 

3.67 ±0.12 

4.33 ± 0.05 3.2+0-S 

4.89 ± 0.08 

T 

9.5 ± 0.2 

6.0 ± 0.2 

9.3 ±0.2 

8.98 ± 0.09 9.6+°| 

6.0 ±0.2 

NcompTT (10 2 > 


6 7 +0 - 5 
- 0.4 



6.7 ±0.6 


Xi/d.o.f.) 


1.12(2493) 



1.04(2463) 





Model: cyclabs*Comptt 



N H (10 22 cm- 2 ) 


0.100 ±0.006 



0.104 ±0.004 


Nh pc (10 22 cm- 2 ) 


4.54 ± 0.08 



4 tc+0.16 


/ 


0.645 ± 0.005 



0.643 ± 0.005 


EEdgel (keV) 


7 20 +0-02 

' '^- 0.04 



7 20 +0 02 


TEdgel 

0.075 ± 0.013 

0.125 ± 0.011 0.075 ±0.013 

< 0.084 

0.069 ±0.013 

0.122 ±0.011 0.069 ±0.013 

< 0.084 

EEdge2 (keV) 


8.43 ± 0.05 



8.43 ± 0.05 


TEdge 2 

0.10 ±0.02 

0.078 ±0.013 0.10 ±0.02 

0 i2 +0 - 08 
U ' 1Z, -0.10 

0.10 ±0.02 

0.079 ±0.013 0.10 ±0.02 

0 i2 +0 08 

U,A -0.10 

Ecyclabs (keV) 


0 69 +0 03 
- 0.04 


0.68^ 

W C y C iabs (keV) 


0.16 ±0.04 



0.15 ±0.04 


Deptcydabs 


0.13 ±0.02 



0.12 ±0.02 


kT 0 (keV) 


0.055 ± 0.010 



0.051 ±0.011 


kT e (keV) 

3.64 ±0.13 

4.34 ± 0.05 3.2)°-5 

8.91 ±0.10 9.5+° : | 

4.89 ± 0.08 

3 70 +0 - 13 
' -0.11 

9.2 ±0.2 

4.34 ± 0.05 3.2 + g! 

4.89 ± 0.08 

T 

9.4 ±0.2 

6.0 ± 0.2 

8 9i+ 0 -io 9 ^+ 0.4 

6.0 ±0.2 

NcompTT (. 1U ^ J 


6.8 ±0.5 



6.9 ±0.6 


xl/ d -°-f) 


1.12(2493) 



1.04(2463) 



Note —Uncertainties are at the 90% confidence level for a single parameter. Th e xl ed (d.°.f.) values are obtained taking the emission lines shown 
in Table[3]into account. 


by llaria et ah (1201 31. Finally, we added two absorption edges 
at 7.2 and 8.4 keV. Fitting the spectra we obtain a^f 2 (d.o.f.) of 
3024(2496) and large residuals between 0.35 and 1 keV are vis¬ 
ible. We show the residuals in Fig.[8](top panel), the best-fit val¬ 
ues of the continuum emission and absorption edges in Tabled 
and the best-fit parameters associated with the emission lines in 

Tabled 


To fit the large residuals between 0.4 and 1 keV, we added 
a Gaussian absorption line (gabs in XSPEC) that is a multi¬ 
plicative component. The component gabs is defined by three 
parameters, and its functional form is 

M(E) - exp[-(r/ V27r<x) exp(-((£ - E 0 )/4cr)) 2 ], 

where r, cr, and Eq are the line depth, the line width in keV, 
and the line energy in keV, respectively. We interpret this 
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Table 3. Best-fit values of the emission lines 



MOS12 spectrum in 0.6-7 keV 

MOS12 spectrum in 1.5-7 keV 

Line 

E (keV) 

o-(eV) 

I (xl0~ 4 ) 

E (keV) 

cr (eV) 

I (xlO -4 ) 




Model: 

Comptt 



N vii 

0.5 (fixed) 

1 (fixed) 

2.8 ± 0.9 

0.5 (fixed) 

1 (fixed) 

2.5 ± 0.9 

O vii (i) 

0.5687 (fixed) 

1.5 (fixed) 

19 ±2 

0.5687 (fixed) 

1.5 (fixed) 

17 ±2 

O viii 

0.6536 (fixed) 

1 (fixed) 

1.8 ±0.5 

0.6536 (fixed) 

1 (fixed) 

1.6 ±0.5 

Ne ix (i) 

0.915 (fixed) 

3 (fixed) 

2.5 ± 0.3 

0.915 (fixed) 

3(fixed) 

2 6 +0 - 9 
-Q-3 

Nex 

1.022 (fixed) 

3 (fixed) 

0.9 ±0.2 

1.022 (fixed) 

3(fixed) 

0 9+0- 5 
-0.2 

Mg xi (i) 

1.3434 (fixed) 

3 (fixed) 

0.72 ±0.14 

1.3434 (fixed) 

3(fixed) 

0 69+°* 8 
-Q-13 

Mg xii 

1.4726 (fixed) 

3 (fixed) 

0.40 ± 0.06 

1.4726 (fixed) 

3(fixed) 

0 34+°-^ 
-Q-Q2 

Si xiv 

2.005 (fixed) 

3 (fixed) 

0.35 ±0.10 

2.005 (fixed) 

3(fixed) 

0 35+0-2 7 

J -0.08 

Fei 

6.408 ± 0.005 

20 ± 15 

2.31 ±0.15 

6.408 ± 0.005 

22 ± 15 

2.30 ±0.15 

Fe xxvi 

6.98 ± 0.03 

50 (fixed) 

0.59 ±0.13 

6.98 ± 0.03 

50 (fixed) 

0.60 ±0.14 




Model: gabs*Comptt 



N vii 

0.5 (fixed) 

1 (fixed) 

1.7 ±0.8 

0.5 (fixed) 

1 (fixed) 

1.7 ±0.8 

Ovu (i) 

0.5687 (fixed) 

1.5 (fixed) 

14.8 ± 1.5 

0.5687 (fixed) 

1.5 (fixed) 

14.6 ±1.5 

O viii 

0.6536 (fixed) 

1 (fixed) 

1.8 ±0.5 

0.6536 (fixed) 

1 (fixed) 

1.8 ±0.5 

Neix (i) 

0.915 (fixed) 

3 (fixed) 

2.6 ± 0.4 

0.915 (fixed) 

3(fixed) 

2 - 6 -m 

Nex 

Mg xi (i) 

1.022 (fixed) 
1.3434 (fixed) 

3 (fixed) 

3 (fixed) 

0.7 ± 0.2 

0.68 ±0.14 

1.022 (fixed) 
1.3434 (fixed) 

3(fixed) 

3(fixed) 

0 66 + ^ 56 
-Q-13 

Mg xii 

1.4726 (fixed) 

3 (fixed) 

0.39 ± 0.06 

1.4726 (fixed) 

3(fixed) 

0 34+0-21 

-Q1Q 

Si xiv 

2.005 (fixed) 

3 (fixed) 

0.40 ± 0.05 

2.005 (fixed) 

3(fixed) 

0 36+° ^ 
^•^^-0.08 

Fei 

6.408 ± 0.003 

23 ± 10 

2.4 ± 0.2 

6.408 ± 0.005 

24 ± 15 

2.3 ±0.2 

Fe xxvi 

6.98 ± 0.02 

50 (fixed) 

0.62 ±0.13 

6.98 ± 0.03 

50 (fixed) 

0.62 ±0.14 




Model: cyclabs*Comptt 



N vii 

0.5 (fixed) 

1 (fixed) 

1.6 ±0.8 

0.5 (fixed) 

1 (fixed) 

1.6 ±0.8 

O vii (i) 

0.5687 (fixed) 

1.5 (fixed) 

14.8 ± 1.5 

0.5687 (fixed) 

1.5 (fixed) 

14.6 ±1.5 

O viii 

0.6536 (fixed) 

1 (fixed) 

1.8 ±0.5 

0.6536 (fixed) 

1 (fixed) 

1.8 ±0.5 

Neix (i) 

0.915 (fixed) 

3 (fixed) 

2.6 ± 0.4 

0.915 (fixed) 

3(fixed) 

2 6 +lo 
a 8 f 

0 63+°-* 7 
-Q-13 

Nex 

1.022 (fixed) 

3 (fixed) 

0.8 ±0.2 

1.022 (fixed) 

3(fixed) 

Mg xi (i) 

1.3434 (fixed) 

3 (fixed) 

0.64 ±0.14 

1.3434 (fixed) 

3(fixed) 

Mg xii 

1.4726 (fixed) 

3 (fixed) 

0.32 ±0.10 

1.4726 (fixed) 

3(fixed) 

0 32+° ^ 
-Q1Q 

Si xiv 

2.005 (fixed) 

3 (fixed) 

0.37 ± 0.09 

2.005 (fixed) 

3(fixed) 

0 37+0-2/ 

Fei 

6.408 ± 0.005 

26 ± 15 

2.4 ± 0.2 

6.408 ± 0.005 

26 ± 15 

2.4 ± 0.2 

Fe xxvi 

6.98 ± 0.03 

50 (fixed) 

0.64 ±0.14 

6.98 ± 0.03 

50 (fixed) 

0.65 ±0.14 


Note —Uncertainties are at the 90% confidence level for a single parameter. The line intensities are in units of photons cm 2 s *. 


component as a CRSF in the spectrum. The model becomes 
Ed*phabs*(f*cabs*phabs*(LN+gC)+(l-f)*(LN+gC)), 
where gC is gabs*CompTT and Ed takes the absorption edges 
into account. The addition of the gabs component improves 
the fit. We obtain a * 2 (d.o.f.) of 2798(2493) with a A\ 2 of 226 
and a F-statistics value of 67.1. The residuals are shown in Fig. 
[8] (the second panel from the top). The best-fit parameters are 
shown in Tables [2] and [3 

We also fitted the residuals between 0.4 and 1 keV using, 
instead of a Gaussian absorpti on line, an absorption line with a 
Lorentzian shape (Mihara et al.lll990l) . This multiplicative com¬ 
ponent (eyelabs in XSPEC) is defined by three parameters, and 
its functional form is 

M(E) = exp[-r (crE/E 0 ) 2 /[(E - E 0 f + cr 2 ]], 

where r, cr, and Eq are the line depth, the line width in keV, and 
the line energy in keV, respectively. We fixed the second har¬ 
monic depth to zero in the model. In this case the adopted model 
is Ed*phabs*(f*cabs*phabs*(LN+cC)+(l-f)*(LN+cC)), 
where cC is cyclabs-'CompTT. The addition of the cyclabs 
component instead of the gabs component gives an equivalent 
fit with a;y 2 (d.o.f.) of 2795(2493). The best-fit values are shown 
in Tables [2] and [3] The residuals are shown in Fig. [3 (bottom 
panel). 


The residuals in Fig.[8]show that the RGS12, XIS1, XIS023, 
MEG, and EPIC-pn spectra are in good agreement below 1.5 
keV, unlike in the MOS12 spectrum. For this reason we repeated 
the analysis described above excluding the 0.6-1.5 keV energy 
range in the MOS12 spectrum. Fitting the spectra using the ini¬ 
tial model, we find a / y 2 (d.o.f.) of 2713(2466) and large residu¬ 
als between 0.35 and 1 keV are visible (see Fig. [9] top panel). 
Using the model that includes the gabs component, we obtain 
a;y 2 (d.o.f.) of 2565(2463) with a Ay- 2 of 148 and a F-statistics 
value of 47.4. The residuals are shown in Fig. [9] (middle panel). 
Using the model including the cyclabs component, we obtain 
a ^ 2 (d.o.f.) of 2565(2463). The residuals are shown in Fig. [9] 
(bottom panel). The unfolded spectra relative to the model in¬ 
cluding the gabs component is shown in Fig. HO] We show the 
data/model ratio with respect to the best-fit model, but excluding 
the gabs component in Fig.ITTI 


Furthermore, using the initial model with the MOS12 spec¬ 
trum between 1.5-7 keV, we look for the presence of CRSF at 
high energies, as suggested by ISasano et al.l (i2014l) . exploiting 
the availability of the JEM-X and IBIS spectra. Adding to the 
model a cyclabs component with width and centroid fixed to 5 
keV and 33 keV (see ISasano et al.lf2014l) . respectively, we find 
an upper limit on the depth of 0.10 at 99.7% confidence level 
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Fig. 9. Residuals with respect to the best-fit models shown in 
Tables [2] and [3] The colours are defined as in other figures. The 
data are graphically rebinned. From top to bottom, the residuals 
with respect to the continuum consist of: 1) a Comptt partially 
absorbed by local neutral matter (large residuals are evident at 
0.7 keV); 2) gabs*Comptt with the energy of the gabs compo¬ 
nent close to 0.72 keV; 3) cyclabs*Comptt. The MOS12 spec¬ 
trum covers the 1.5-7 keV energy band. 

(3<x). This val ue i s not consistent with the 0.4^ obtained by 
ISasano et alJ (120 14l> . 

We find that the energy and width of the CRSF at 0.7 keV are 
consistent for the Lorentzian and the Gaussian shapes. The depth 
values of the absorption edges at 7.2 keV are compatible for the 
XMM-Newton and INTEGRAL spectra, while they higher in the 
Suzaku spectrum. Finally, the optical depth of the Comptonised 
component assumes similar values close to nine in the XMM- 
Newton, Chandra, and Suzaku spectra, while it is 6.0 + 0.2 in the 
INTEGRAL spectrum. 

Finally, we note that including MOS12 spectrum between 
0.6 and 1.5 keV only marginally affects the fit results. In the 
following we use the best-fit values obtained when excluding 
the 0.6-1.5 keV energy range of MOS12. 

6. Discussion 

We used three non-simultaneous pointed X-ray observations of 
X1822-371: a XMM-Newton observation (using RGS, MOS, 
and EPIC-pn spectra), a Suzaku observation (using XIS and 
HXD/PIN spectra), and finally, a Chandra observation (using 
the first-order MEG spectrum). The spectral analysis of the 
XMM and Chandra data sets have already been discussed in 
llaria et a i](l2013l). The same XMM/EPIC-pn data were also anal¬ 
ysed IryTsomenTetak] d2012l) . who report a different interpre¬ 
tation of the X-ray spectrum. Moreover, we used all the avail¬ 
able INTEGRAL/JEM-X and INTEGRAL/ISGRI observations 
of XI822-371 and extracted the corresponding spectra to con¬ 
firm or disprove a c laim ed CRSF at 33 keV in the Suzaku/PIN 
spectrum (see lSasano et al.ll2014l) . 


Fig. 10. Unfolded spectra relative to the model including the 
gabs component. The MOS 12 spectrum ranges between 1.5 and 
7 keV. Colours as above. 



Fig. 11. Data/model ratio with respect to the best-fit model 
shown in Tables[2]and[3] but excluding the gabs component. The 
data are graphically rebinned. Colours as above. The MOS 12 
spectrum ranges between 1.5 and 7 keV. 


We adopted the same model as proposed by those authors to 
fit the continuum emission. It consists of a Comptonised com¬ 
ponent Comptt absorbed by interstellar matter and partially ab¬ 
sorbed by local neutral matter. The effect of Thomson scatter¬ 
ing on the local cold absorber is taken into account by adding 
the cabs component with an equivalent hydrogen column den- 
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sity imposed t o be the sa me as the one associated with the local 
neutral matter. ITaria et al. . ( 2013 3 found residuals in the EPIC-pn 
data between 0.6 and 0.8 keV and added a black-body compo¬ 
nent with temperature fixed to 0.06 keV, thereby improving the 
fit significantly. 


In this work we give a different interpretation of the residuals 
in the EPIC-pn data between 0.4 and 1 keV, modelling them with 
the addition of a CRSF close to 0.7 keV. The broad residuals are 
fitted using the gabs or the cyclabs model components. The 
addition of this component significantly improves the fit, and we 
obtain statistically equivalent fits using either the gabs or the 
cyclabs component. 

The interpretation of the residuals as a CR SF in the spect rum 
allows us to refine the scenario proposed bv ilaria et al.l (12013h for 
XI822-371. The authors suggested that the Comptonised com¬ 
ponent originates in the inner region of the system, it is not di¬ 
rectly observable because of the large inclination angle of the 
system, and only 1 % of its flux arrives to the observer because 
of scattering by an extended optically thin corona with optical 
depth ~ 0.01. We now suggest that the Comptonised component 
could be produced in the accretion column onto the NS magnetic 
caps. 


Below we discuss our results showing that a CRSF at 0.7 
keV agrees with the non-conservative mass transfer scenario 
proposed for XI822-371 in the past three years by several au- 
thors (see, for example, Iaria et all 1201 3l 1201 lb IBurderi et alJ 
1201 Ot IBavless et al.l I201C ) and allows the NS mass of the bi¬ 
nary syste m that is between 1.61 and 2.32 M 0 to be furthe r 
constrained! Munoz-Darias et ap 120051). IBurderi et al.l j20lol) . 
IBavless et ah ( 201 Oil , and Iaria et al. ( 201 lb found a large orbital 
period derivative of XI 822-371, and Burderi et al.l (1201 Oh show 
that this indicates that XI822-371 accretes at the Eddington limit 
and that the rest of the mass transferred by the companion star is 
expelled from the system. 


Adopting the values for P s and P s derived in Sect. [3] and as¬ 
suming that XI822-371 accretes at its Eddington limit, we show 
that the CRSF energy obtained by our fits is consistent with a 
scenario in which the NS in X1822-371 is spinning up. Since 
P s is negative, the corotation radius r c , the radius at which the 
accretion disc has the same angular velocity as the NS, has to be 
larger than the magnetospheric radius r m , the radius at which the 
magnetic pressure of the NS B-field equals the ram pressure of 
the accreting matter. The corotation radius can be expressed as 


r c = (GM^s , where M N $ is the NS mass and P s is 

the NS spin period. Assuming a NS mass of 1.4 M 0 , we obtain 
r c - 1,180 km; for a NS mass of 2 M Q , the corotation radius is 
r c =: 1,300 km. The magnetospheric radius is given by the re¬ 
lation r m = (hr A , where r A js the Alfven radius and (f> a constant 
close to 0.5 (see Ghosh & Lambll99llf . We estimate r A using the 
Eq. 2 in IBurderi et al.l (1 19981) : 


r A = 4.3 x Rf n L~ 2/1 e 2n m' n km, (1) 

where //30 is the magnetic moment in units of 10 30 G cm 3 , Re the 
NS radius in units of 10 6 cm, L 37 the luminosity of the system in 
units of 10 37 erg/s, m the NS mass in units of solar masses, and 
finally, e is the ratio between the luminosity and the total grav¬ 
itational potential energy released per second by the accreting 
matter. 

Adopting the best-fit value of the CRSF energy, E^v, ob¬ 
tained from the gabs component, and assuming that this is pro¬ 
duced at (or very close to) the NS surface, we estimate the NS 


B-field using the relation Ek ev = (1 + z) 1 11.6/1 12 , where 


d+zr 1 


/ 2 GM ns \ U2 

l Rnsc 2 ) 


and B 12 is the NS B-field in units of 10 12 G. We arrange the term 
(1 + z) _I in terms of m and R(, and obtain 

fl, 2 = ts( i_o295 £) ,,2 °' (2 > 

We assume that the intrinsic luminosity of X1822-371 is the 
Eddington luminosity, L - 1.26 x 10 38 (M ra /M 0 ) erg/s, which 
we rewrite as L 37 = 12.6 m, which is the Eddington luminosity 
in units of 10 37 erg/s. Substituting this expression of luminosity 
and the expression of Bn of Eq.[2]into Eq.|T|and assuming a NS 
radius of 10 km, we obtain 

r A = 4.3 x 10 2 (1 - 0.295m)“ 2/7 »r 1/7 km. 


For NS masses of 1.4 and 2 M e , we obtain r A - 490 and ^510 
km, respectively. Since r m - cf>r A with </> - 0.5 then for a NS 
mass of 1.4 and 2 M 0 , we obtain r m - 245 and - 255 km, 
respectively. This means that, for a NS mass between 1.4 and 
2 M 0 , r m is always smaller than r c by a factor five. This implies 
that the accreting matter gives specific angular momentum to the 
NS, which increases its angular velocity and spin s it up. _ 

Ne xt we adopt the set of relations shown bv lGhosh & Lambl 
( 1979 ) (Eqs. from 15 to 18) to establish a relation for the deriva¬ 
tive of the spin period, the spin period, the luminosity, and the 
NS mass. This set of equations is valid for the fastness pa¬ 
rameter a> s — D-JUk^o) < c^>max — 0.95, where Cl s is the 
NS angular velocity, il K (r 0 ) is the Keplerian angular velocity 
at r (> , and rg is the radius that separates the boundary layer 
from the outer transition zone (see discussion in lGhosh & Lambl 
11979h . Using Eq. [2] and assuming that the NS accretes at the 
Eddington limit (L 37 = 12.6 /«), the parameter co s of the Eq. 16 
in lGhosh & Lambl ( 1979 ) becomes 

Wi "°- 456 (ff^) ( 1_0 - 295 ^) 5P >n- 5n P^. (3) 

Using the value of E[ :i , v shown in Table [2] the spin period value 
shown in Sect. [3] and finally, imposing that Rg — 1, we find that 
w s is between 0.063 and 0.083 for m between 1 and 3 M 0 . This 
implies that the values of NS B-field and luminosity satisfy the 
spin-up condition. 

Using Eq. 15 of Ghosh & Lambl (1 19791) . we constrained the 
NS mass . We adopted the expressio n of the NS moment of inertia 
given bv lLattimer & Schutzl (120051) in Eq. 16 of their work. The 
expression is valid for many NS equation of states and for a NS 
mass higher than 1 M 0 . Rewriting the expression in terms of m 
and R(, we obtain 


I 45 * (0.471 ± 0.016) mRl 


1 +0.42— +0.009^-1, 

Re Rt! 


(4) 


where 745 is th e NS mom ent of inertia in units of 10 45 g cm 2 . The 
Eq. 15 of lGhosh & Lambl ( fl979l) can be rewritten as 


- P 12 =* 29.51 C 2/7 m- 4/7 ^V 1 / ,2 n(w J ), 


(5) 


where P -12 is the spin period derivative in units of 10 12 s/s, 6 
is the term in parenthesis in Eq. [4] B n is given by Eq. [2] and 
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finally, the function n(u> s ) in its useful approximate expression 
is 


n{oj s ) ~ 1.39{1 - w,[4.03(l - w,) 0173 - 0.878]}(1 - w,)' 1 

(see Eq. 10 in lGhosh & Lamblll979l) . 

Initially, we impose that R(, = 1 and find how P. 12 changes 
as function of m. The error associated with P. 12 is mainly due to 



Fig. 12. Spin period derivative, P- 12 , vs. m for a NS radius 
of 10 km (red curve); the black curves represent the upper 
and lower limit values of the spin derivative. The grey box in- 
dicates the mass range al lowed for XI822-371 according to 
iMunoz-Darias et akl ( i2005l) . The horizontal blue strip indicates 
the value of P -12 as derived in Sect. [3] 


the term n(a> s ), that has an accuracy of 5% fsee lGhosh & Lambl 
fl979h : we take also the errors associated with I and Ek e v into 
account. We adopt the lcr error for Ek e v- We show the values of 
P 12 vs in in Fig. |T2] For a range of m between 1.3 and 2.5, the 
P -12 changes from -3.4 up to -1.6. The grey box in Fig.fTSllimits 
the allowed NS mass betw een 1.61 and 2.32 M Q , according to 
IMunoz-Darias et afl (2005 1. Furthermore, the value of the spin 
period derivative, P-n = -2.55 + 0.03 s/s that we obtained in 
Sect. [3] is shown, as are the uncertainties associated with P-12. 
We find a NS mass of 1.69 + 0.13 M n assuming^ = 1 , whic h 
is inside the range suggested by [Munoz-Darias et akl (20051) . 
Consequently, using the Eq. [2j we find that the NS B-field is 
(8.8 + 0.3 ) x IQ 10 G , a v alue that is very similar to the one sug¬ 
gested by iJonker & van der Klisl (200 ll) assuming a luminosity 
of ~ 10 38 erg/s for XI822-371. 

Using the mass function (2.03 ± 0.03) X 10 ~ 2 M 0 and an in¬ 
clination angle of X1822-371 of 8235 fsee IJonker et alJ [20031) . 
we infer the companion star mass, M c - 0.46 + 0.0 2 M f . ; , which 
is clos e to the value of 0.5 M 0 suggeste d by Munoz-Darias et akl 
(20051) . Using only optical observations, [SomeroetaL ( 2012 ) es¬ 
timate that the mass ratio of X1822-371 is q — M r j M N s = 0.28. 
Using this relation, we find that the companion star mass is 
M c - 0.47 ± 0.04 M 0 for a NS mass of 1.69 + 0.13 M Q , and 


that value is compatible with the one inferred by us using the 
mass function. 


Table 4. Values of IUs, Myy, B, and M c obtained from the cy¬ 
clotron line energy found adopting the gabs component in the 
fit. 


Rms 

(km) 

Mjvs 

(M 0 ) 

B 

(10 10 G) 

M, 

(M 0 ) 

8 

1.70 + 0.11 

10.2 + 0.3 

0.46 ± 0.02 

8.5 

1.71+0.12 

9.7 ± 0.3 

0.46 ± 0.02 

9 

1.71+0.12 

9.4 ± 0.3 

0.46 ± 0.02 

9.5 

1.70 + 0.13 

9.0 ±0.3 

0.46 ± 0.02 

10 

1.69 + 0.13 

8.8 ±0.3 

0.46 ± 0.02 

10.5 

1.67 + 0.14 

8.5 ± 0.3 

0.45 ± 0.03 

11 

1.64 + 0.14 

8.3 ±0.3 

0.45 ± 0.03 

11.5 

1.61+0.15 

8.1 ±0.3 

0.44 ± 0.03 


Note -Uncertainties are discussed in the text. For clarity we also show 
the values of My,;, B, and M (; for R N4 - = 10 km. The errors are at 68% 
confidence level. 


Assuming different NS radii for Eq. [5] we obtain different 
values of the NS mass. We show the NS masses for several val¬ 
ues of the NS radius ranging from 8 to 11.5 km in Table[4] The 
NS radius in X1822-371 cannot be larger than 11.5 km because 
the NS mass w ould be lower than 1. 61 M 0 which is the lower 
limit given by IMunoz-Darias et all ( 20051) . This result further 
constrains the NS mass range between 1.46 and 1.81 M 0 . The 
NS B-field ranges between 7.8 and 10.5 xlO 10 G, and finally, 
the companion star mass ranges between 0.41 and 0.48 M 0 . 

We comp ared the results in Table [4] with those obtained by 
ISteiner et al.l (20101) , which determined an empirical dense mat¬ 
ter equation of state from a heterogeneous data set of six neutron 
stars: three Type-I X-ray busters with photospheric radius expan¬ 
sion and three transient low-mass X-ray binaries. O ur compari- 
son was done with the results reported in Table 7 bv ISteiner et al.l 
( 2010 1). They are valid for a NS radius equal to the photospheric 
radius of the NS. The authors find that for a NS mass of 1.6, 
1.7, and 1.8 M e , the corresponding radius is 10.8^®, 10.7/j 8 , 
and 10.7 /j j km, respectively, with the errors at 95%. If the 
NS i n X1822 -371 is similar to those of the sample studied by 
ISteiner et al.l ( 20101) . we can exclude from Table |3] the solutions 
for NS radii smaller than 9.5 km. This implies that the NS mass 
range is between 1.61 + 0.15 and 1.70 ± 0.13 M Q , that the NS 
B-ffeld in units of 10 10 G is between 8.1 + 0.3 and 9.0 + 0.3, and 
finally that the companion star mass is between 0.44 ± 0.03 and 
0.46 ± 0.02 Mo. 

We note that the estimation of the CRSF energy is model de¬ 
pendent. The CRSF energy is 0.72 and 0.68 keV, adopting the 
gabs and cyclabs components, respectively, to fit the averaged 
spectrum. However, the values of the NS mass, NS magnetic- 
field strength, and companion star mass at different NS radii are 
the same as shown in Table [4] even using the CRSF energy ob¬ 
tained from the cyclabs component, since this energy and the 
NS magnetic-field strength, B\ 2 , are linearly dependent (see Eq. 
13 and because the spin period derivative weakly depends on B \2 
(see Eq.|3- 

The gabs component allows us to estimate the temperature 
of the plasma where the CRSF originates, assuming that the 
broadening of the line has a thermal origin. At the cyclotron res¬ 
onance frequency a> c , electrons at rest absorb photons of energy 
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Hu> r . Fo r thermal Doppler broadening, Aoj d is predicted to be 
dMeszaroJl992h 


A cod 
oj c 



| cos 6\, 


where hAajo - cr ga b s , hoj c = Eg ahs , kT is the electron temper¬ 
ature, and m e c 2 is the electron rest energy. The angle 6 mea¬ 
sures the direction of the magnetic field with respect to the line 
of sight. Outside the range a> c ± A a>o, the cyclotron absorption 
coefficient decays exponentially, and other radiative processes 
become important. Substituting the values of cr ga b s and £ ga b s 
shown in Table [2] (for MOS12 spectrum ranging between 1.5 
and 10 keV), we obtain a lower limit on the plasma temperature 
of kT = 8 + 2 keV with the error at 68% confidence level, which 
is a factor of two or three larger than the electron temperature of 
the Comptt component, but the values are consistent at the 2cr 
level. This suggests that the Comptonised component is probably 
produced in the accretion column onto the NS magnetic caps. 

We do not observe cyclotron harmonics in the spectrum. 
To date, the lowest energy measured for a CRSF produced by 
electron motion around the NS magnetic;field line s is 9 keV in 
the source XMMU J054134.7-682550 dManousakis et alj[2009h : 
also in that case, harmonics are not visible in the spectrum. The 
cyclotron line observed in the spectrum of the Be/X-Ray Binary 
Swift J1626.6-5156 has an energy of 10 keV and only_a weak in¬ 
dication of a harmonic at 19 keV (see lDeCesar et al.ll2013lf . The 
source KS 1947+300 has been recently observed with Nuclear 
Spectroscopy Telescope Array ( NuSTAR ) and Swift/XRT in the 
0.8-79 keV energy range (tFiirst et al.il2014h : a CRSF at 12.5 keV 
has been observed but no harmonics have been detected. Finally, 
the anomalous X-ray pulsar SGR 0418+5729 shows a CRSF 
produced by proton motion,_with centroid at 1 keV, and no har¬ 
monics are observed (see lTiengo et al.ll2013 ij). To now, only the 
isolated NS CCO 1E1207.4-5209 shows a CRSF also and its first 
ha rmon i c at 0.7 and 1.4 keV, respectively dSanwal et al.l 12002 : 
iMereghetti et al. 2002). These results show that, although pecu¬ 
liar, it is possible to observe only the fundamental harmonic of 
the CRSF. 

Finally, we note that the presence of a CRSF at 0.7 keV 
shown in this work contrasts with the recent result suggested by 
ISasano et al.l t 2014 ). who find a CRSF at 33 keV when analysing 
the same Suzaku data as presented in this work. First of all, we 
note that a CRSF at 33 keV is detectable including the HXD/PIN 
data up to 40 keV. However, we have shown that the HXD/PIN 
source spectrum is overwhelmed by the NXB+CXB spectrum at 
energies higher than 36 keV (see Fig. [7j. Furthermore, we note 
that, at 33 keV, the HXD/PIN effective area is nearly 50 cm 2 , 
while the HPGSPC and PDS instruments on-board BeppoSAX 
had an effective area of ~ 200 and ~ 500 cm 2 , respectively. 
Ilaria et ahld200llf analysed a broad band spectrum of XI822-371 
using the narrow-field instruments on-board BeppoSAX and did 
not find any evidence of a CRSF at 30 keV with a PDS exposure 
time of 18.7 ks. The HXD/PIN spectrum has an exposure time 
of 37.7 ks, which is a factor of two longer than the PDS expo¬ 
sure times, but it has an effective area a factor of 10 smaller at 
33 keV. The presence of a CRSF at 33 keV in the Suzaku data 
is therefore unrealistic when assuming that it does not change in 
time. 

To verify the presence of a CRSF at 33 keV, we have also 
analysed the IBIS and JEM-X spectra for an effective dead-time- 
corrected exposure of 874 ks and 283 ks, respectively. We fit¬ 
ted both the spectra, adopting the same model as used to fit the 


XMM-Newton/Suzaku data and found no evidence of CRSF at 
33 keV. 

We also note that a CRSF at 33 keV is not consistent with the 
observed spin-up of the NS. Assuming that the NS is spinning- 
up and using the NS B-field value inferred by a CRSF at 33 
keV, ISasano et al.l ( 2014 1 obtain an intrinsic luminosity of the 
system of ~ 3 x 10 37 erg s _1 needed to have the measured spin- 
up rate. In case of spin-up, we expect that the corotation radius 
r c has to be larger than the magnetospheric radius r m . The values 
of r c is 1 180 km and 1 300 km for a NS mass of 1.4 and 2 
M 0 , re spectively, assuming the value reported by ISasano et al.l 
( 201411 as spin period. To estimate the magnetospheric radius, 
we used the Eq. Q] assuming a NS radius of 10 km and a NS 
B-field of 2.8 x 10 12 G (see Sasano et al.ll2014l) . For a NS mass 
of 1.4 M 0 , we find that r A - 5 940 km and r m - 4>r A - 3 000 
km, with 4> - 0.5. We infer that r m = 3r c , and this implies that 
matter cannot acc rete onto the NS ; this s cenario corresponds to 
that discussed bv iGhosh & Fambl (1979) for ai s » (o max with 
(Omax - 0.95. In fact, 

co s = l.35B 6 1 ,2 m- 2ll p- l L^ n 

for a NS radius of 10 km (see Eqs. 16 and 18 in lGhosh & Fambl 
Il979lf . Using th e values of luminosi ty, spin period, and NS B- 
field reported bv ISasano et al] (l2014lf . we find that a> s — 3.2 for 
a NS mass of 1.3 M 0 and = 2.5 for 3 M e . This result suggests 
that the scenario is not self-consistent because the values of lu¬ 
minosity and NS B-field would contradict the observed spin-up 
of the NS. 


7. Conclusion 

We analysed the broadband X-ray spectrum of X1822-371 using 
all the presently available data sets that allow for high-resolution 
spectroscopic studies. Our aim was to understand the nature of 
the residuals between 0.6 and 0.8 keV previously observed in the 
XMM/EPIC-pn data bv llaria et al. ( 2013 ). Initially, we analysed 
the averaged spectra simultaneously. The broad band spectrum 
obtained from Suzaku constrains the continuum emission well. 
The adopted model was a Comptonised component absorbed 
by the neutral interstellar matter and partially absorbed by local 
neutral matter. We took the Thomson scattering into account and 
fixed the value of the equivalent hydrogen column density of the 
local neutral matter to that producing the Thomson scattering. 

The values of the parameters are similar for the three data 
sets and consistent with previous studies (see e.g. Ilaria et akl 
120131) . To fit the residuals between 0.6 and 0.8 keV, we added 
an absorption feature with Gaussian profile (gabs in XSPEC). 
Alternatively, we adopted an absorbing feature with Forentzian 
profile (cyclabs in XSPEC). In both cases the addition of a 
CRSF to the model improved the fit. We found that the improve¬ 
ment does not depend sensitively on the exact shape used to 
model the absorption profile. 

We also detected the spin period of X1822-371 in the EPIC- 
pn data. We obtained the value of 0.5928850(6) s. Using all 
the measurements known of the spin period of XI822-371, 
we estimated that the spin period derivative of the source is 
-2.55(3) x 10- 12 s/s, and this confirms that the neutron star is 
spinning up. Folding the EPIC-pn light curve, we derived a pulse 
fraction of 0.75% in the 2-5.4 keV energy band. 

Using the best-fit values of the CRSF parameters, under the 
assumption that the system is accreting at the Eddington limit, 
we estimate a NS B-field between (8.1 + 0.3) x 10 10 and (9.0 + 
0.3) x 10 10 G for a NS radius ranging between 9.5 and 11.5 km. 
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We subsequently constrain the NS mass assuming that the CRSF 
is produced at the NS surface. We find that, for a Gaussian profile 
of the CRSF, the NS mass is between 1.61 ±0.15 and 1.70 + 0.13 
M 0 . The companion star mass is constrained between 0.44+0.03 
and 0.46 ± 0.02 M Q . 

Finally, we note thatour conclusions contrast with the recent 
results reported by ISasano et al.l (120141) . who report detecting a 
CRSF at 33 keV (and a corresponding NS-B field of 3 xlO 12 
G) in the Suzaku data also used in this work. To address this 
point, we have selected the whole IBIS and JEM-X public data 
set of the XI822-371 region. We extracted the JEM-X and IBIS 
spectra of X1822-371 having an exposure time of 330 and 874 
ks, respectively. The INTEGRAL spectrum combined with the 
XMM, Suzaku, and Chandra spectra does not show a CRSF at 
33 keV. Furthermore, we also show from theoretical arguments 
that a CRSF at 33 keV is not consistent with the evidence that 
the NS in X1822-371 is spinning up. 
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